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Abstract 

Ridge-like correlations in high-energy proton-proton collisions reported by the 
CMS collaboration suggest a collective flow that resembles the one in heavy-ion 
collisions. If the hydrodynamic description is valid then the effect results from the 
initial anisotropy of the colliding matter which depends on the structure of protons. 
Following recent theoretical developments, we propose several phenomenological 
models of the proton structure and calculate the anisotropy coefficients using the 
Monte Carlo Glauber model. Our estimates suggest that the event multiplicity 
dependence allows one to discriminate between different proton models. 
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1 Introduction 


The analysis of two-particle angular correlations in pp collisions at -y/s = 7 TeV revealed 
an unexpected near side (A 0 = 0 ) correlation in the azimuthal angle of produced particles 
(’ridge effect’) [1]. There is no obvious reason why such a long-range in pseudorapidity 
correlation should occur [2,3]. Ridge effect has been previously observed in relativistic 
heavy-ion collisions and it was explained by a collective anisotropic flow of hot and dense 
medium (hadronic gas or quark gluon plasma) created at the collision. 

It is possible that the hydrodynamic explanation applies to high-multiplicity pp events 
which are sufficiently energetic for the collision products to thermalize [4]. In this paper 
we assume validity of such a scenario and postulate that the proton internal structure 
can manifest itself in the ridge effect. During the hydrodynamic expansion, the spatial 
anisotropies of the collision area are transformed into anisotropic flow of produced par¬ 
ticles. Thus, the spatial anisotropies due to the proton structure should be reflected in 
the measurable collective flow coefficients. 

Actually, there exist in literature estimates of the possible elliptic flow in pp colli¬ 
sions [5-13]. We basically follow the approach of [5] in order to make a link between the 
proton models considered by us and the experimental data. Within the framework of 
Glauber model [14], we calculate for each model the initial spatial anisotropies: eccen¬ 
tricity, denoted by 62 , and triangularity, denoted by 63 , in proton-proton collision taking 
into account event-by-event fluctuations in the proton density prohle. 

The considered phenomenological proton densities are inspired by the quark-diquark [17] 
and RGPEP [18] models of proton. We find that it is possible to distinguish highly ec¬ 
centric (rod-like) and highly triangular (triangle-like) hxed proton densities by looking at 
distributions of 62 and 63 and the multiplicity dependence of their RMS values. We also 
calculate these quantities for the fluctuating anisotropic proton configurations (Gaussian- 
fluctuating model). 

We do not explicitly present results for the anisotropic flow coefficients because 
we limit our discussion to comparison of predictions of different models for the initial 
states in hydrodynamic evolution and the prediction of the exact shapes of the resulting 
final flow anisotropies is beyond the scope of this work. Our calculations concern the 
initial spatial anisotropies, though we emphasize that they should correspond to the flow 
coefficients measurable in multi-particle correlations in a fashion that is unlikely to be 
far from linear [15,16]. 

2 Glauber model for pp collisions 

The optical limit of Glauber model treats the collision of two composite particles as 
a superposition of independent binary collisions between their constituents. In case of 
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pp collisions, we assume the constituents are partons of one type, and their smoothed 
distribution is given by proton density prohle. We £x the partonic cross section agg = 
4.3 mb [15] and adjust the number of partons in proton Ng so that the experimental 
inelastic pp cross section app ~ 60 mb [19] is reproduced. 

The geometrical quantities calculated within the Glauber formalism correspond to 
the impact plane of the collision {xy) such that any ^-dependence is integrated out. The 
density of binary collisions is given by the expression: 


f b 

^co 1 i( 2 :, 2 /;&,Sa, Ss) = agg J dz p \x - 




dz' p[ X +-,y,z'-,'EB ) , ( 1 ) 


where p{x, y, z; S) is a proton density prohle, b is an impact parameter and by S we mean 
a set of huctuating parameters describing the proton internal structure and orientation 
in space (for hxed conhguration models S dependence will correspond only to possible 
rotations in space) [14]. 

For each event characterized by b, and we calculate eccentricity 62 and trian¬ 
gularity 63 . We choose the following dehnition of spatial anisotropy [20]: 


^r). — 
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( 2 ) 


where (j) is an azimuthal angle in xy plane and = x‘^ + y"^. The curly brackets {...} 
stand for the average taken with respect to the impact plane binary collision density 


rr( ^ J dxdy f{x,y)nco\i{x,y]b,i:A,'SB) 
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( 3 ) 


The differential cross section (probability) of a given event depends only on the inte¬ 
grated number of binary collisions Nco\i{b, YIa, ^b) = f dxdyncoii(x,y; b, S^), 


o'ib, Syi, Hb) — 1 
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( 4 ) 


The expectation value of a quantity Q is calculated in the following way 

{Q) = — f ‘^^bdb [ P{'EA)d'EA [ P(SB)dSBcr( 6 ,S^,SB)Q( 6 ,S 2 i,SB), (5) 
J 0 J J 


where app = J^27ibdb J P{^a) dTiA J dTlB o'{b,'SA,'^B) and P(S) is the proba¬ 

bility density of proton conhguration S. 

We calculate charged hadron multiplicity N for each event assuming its linear scaling 
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with A^coii: 


( 6 ) 


N = aA^coii, 

where a is determined by the demand of reproducing experimental minimum bias charged 
hadron multiplicity (N) = 30 [21], 

There is another approach to estimating the multiplicities and geometrical quantities 
within the Glauber formalism, namely the wounded nucleon model [22], Within this 
model, all the averages denoted by are taken with respect to the local density of 

participating constituents instead of the density of binary collisions. However, in high 
energy hadronic collisions, multiple gluonic interactions are expected to occur, so we 
think it is more natural to use here the density of binary partonic collisions. 


3 Models of proton internal structure 


The renormalization group procedure for effective particles (RGPEP) [18] suggests that 
a proton can be described in terms of three effective quarks and a gluon body interacting 
via harmonic potential. For every considered special case of this picture we model the 
proton as a sum of isotropic 3D Gaussian densities pq corresponding to effective quarks 
and anisotropic Gaussians Pg representing the gluonic flux tubes connecting the quarks: 




Pg(r;r„n) 


(2vr)3/2^3® 

1 

(27r)3/2 rj ri 


(7) 

( 8 ) 


The effective quarks and the gluon body are clusters of partons, as in the two stage and 
cluster models [23-25] used to explain the shape of deep inelastic scattering structure 
functions. 


3.1 Fixed configuration models: I and Y 

According to [17], the ridge effect could be the consequence of large eccentricities 62 in 
high-multiplicity collisions between aligned quark-diquark states of protons which con¬ 
stitute around 20 % of all possible states. 

We model this state, labelling it with the symbol I, as a superposition of two effective 
quark bodies, with quark and diquark parton densities generically described in terms 
of simple functions, and two gluonic tubes connecting the quark bodies (Fig. 1). Such 
greatly simplihed density model is sufficient for tracing the effects of symmetric rod-like 
shape of proton in comparison with other conceivable shapes (see below), leaving the 
more subtle effects of the asymmetry between quark and diquark for further analysis 
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elsewhere. Thus, we use here 


pi(r) = N, 
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The free parameters of the model are: the effective quark radius r^, the length of quark- 
diquark tube d, and the percentage of gluon body content n. Assuming that only effective 
quark bodies carry net charge (-|-4/3e homogenously distributed in the diquark and — l/3e 
in the quark), we choose Vq = 0.25fm and d = 1.5fm, which reproduce charge rms radius 
of proton 0.9fm). In this work, we always choose k = 0.5. 

For comparison, we also consider a model that we label by the symbol Y in which 
proton has a highly triangular shape (Fig. 2). The parton density in this case is assumed 
in the form 


PY(r) = Ng 


k=l 


k=l 
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where ri = ^0, ^d, —|j,r 2 = ^0, —^d, —| j , rs = (O, 0, |) and TV(d] is a rotation 
matrix in yz plane. In this model, the effective quarks are located in the vertices of an 
equilateral triangle and the gluon tubes connect them with the center of mass of the 
system. Again we choose Vq = 0.25 fm and d = 1.5 fm. 

The densities of types I and Y are fixed proton structure models in a sense that the 
only parameters that fluctuate event-by-event are the angles describing proton orientation 
relative to the direction of its velocity (not explicitly introduced above). 


3.2 Gaussian-fluctuating conflguration model 

We also consider a model with a fluctuating proton configuration, labelled G-f, in which 
the relative positions of effective quarks differ event-by-event. By generalizing the previ¬ 
ous formulas for proton densities, we introduce Gaussian-fluctuating proton density 

^3 3 

PG-f(r;ri,r2,r3) = Ng^-'^pq{r - rk-,rq) + Ng-'^pg 

k=l k=l 

( 11 ) 

In this expression, TZ[6, 0] transforms vector (0, 0,1) into (cos 0 sin 6, sin 0 sin 6, cos 6) and 
ffc = ?"fc(cos0fc sindfc, sin0fc sin dfc, cosdfc) is the position vector of k-th effective quark. 

Following the harmonic oscillator phenomenology, we assume that the probability 
density of finding a proton whose quarks are in positions ri,r 2 ,r 3 is Gaussian with an 


^[dfc,0fc](r- ^)]rq,^ 
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additional constraint of in the center-of-mass frame (ri + r 2 + ra = 0 , obtained by recen- 
tering of the quarks), 


P(ri,r2,r3) 


^ -(r?+rg+r^)/2fi^p 


( 12 ) 


The free parameters in the G-f model are k, Vg and Rp. Again we set Vg = 0.25 fm, which 
implies the value Rp = 0.43 fm, if one wants to obtain the experimental rms proton 
charge radius. 


4 Results 

We present the event distribution of eccentricities 62 and triangularities 63 for three classes 
of pp collision models: II, YY and G-f in which colliding protons are described by I, Y or 
G-f model, respectively (Fig. 3, 4). We also present the RMS values (e|) and (e|) in 
certain multiplicity bins (Fig. 5, 6 ) and the distribution of event multiplicities N (Fig. 7). 
The reason for presenting RMS and not mean values of anisotropies is because the former 
correspond hydrodynamically to the actual anisotropic flow coefficients extracted from 
two particle correlations (we base this approach on [12]). The minimum bias (averaged 
over all events) \/(^ 2 )mb 0.37, 0.40 and 0.31 for II, YY and G-f model, respectively, 

while a/(^Dmb ~ 0.14, 0.23, 0.15, correspondingly (these quantities can be calculated 
using probability densities from Fig. 3 and Fig. 4). 

All the calculations are performed by means of a Monte Carlo algorithm. 300 000 
events are generated for each class of collision. The statistical weight of an event is given 
by Eq. (4) and the mean values are calculated according to Eq. (5). In case of I and 
Y models, we average only over the impact parameter b and over all the possible space 
orientations of protons during the collision. In case of G-f model, the additional averaging 
over the positions of effective quarks takes place. The event multiplicity N is estimated 
by the use of Eq. ( 6 ). 

We set Ng in such a way that the experimental pp inelastic cross section of 60 mb is 
reproduced. This is accomplished for Ng in range 7-10. Similarly, the parameter a from 
Eq. ( 6 ) is found to vary between 5 and 11. 

5 Conclusions 

Our calculations predict quite high values of eccentricities and triangularities in pp colli¬ 
sions. Assuming the hydrodynamic limit in which the scaling between the initial spatial 
and finite momentum anisotropy is linear and taking ^ 2/62 ~ 0.3 [15], we expect minimum 
bias V 2 ~ 0.11,0.12,0.09 for quark-diquark, triangular and Gaussian-fluctuating model. 
Bozek extracts from the ridge effect V 2 in range of 0.04 - 0.08 [4] which can potentially 
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be in agreement with our results if viscous effects reducing the scaling factor are taken 
into account. 

Our models predict that one in principle should be able to distinguish quark-diquark 
(I), triangular (Y) and Gaussian-fluctuating (G-f) proton conhgurations by extracting 
the collective flow coefficients V 2 and ^3 from multi-particle correlations. The actual 
shapes of the collective flow dependence on multiplicity will differ from Fig. 5 and Fig. 6 , 
though it is reasonable to assume that to a hrst approximation the scaling of Vn with e„ is 
linear. In general, the anisotropies in collisions of triangular protons, denoted as YY, are 
much higher than in those of II or G-f G-f. This was expected because irrespective of the 
orientation of Y protons the collision area is always anisotropic. We hnd that the values 
of anisotropy coefficients in YY collisions increase strongly in high multiplicity bins N. 
Gollisions of quark-diquark protons (II) are characterised by relatively large eccentricities 
62 (though not as large as in YY collisions) and small triangularities 63 . The distributions 
of anisotropies and the N dependence of their mean values are the most smooth for the 
G-f G-f collisions, which results from the existence of additional degrees of freedom. 

We observe that the fraction of high-multiplicity events in the hxed conhguration 
models is smaller than in the Gaussian-fluctuating model. Our comparison to the exper¬ 
imental data (Fig. 7) indicates that a more realistic model of pp collision should include 
not just one type, but all II, YY, G-f G-f (and lY, I G-f, Y G-f) types of collision con¬ 
hgurations with certain probabilities. However, the G-f model alone can approximately 
describe the experimental results for multiplicity distributions. 

Avsar et al. [12] using DIPSY Monte Garlo generator predict values of eccentricity 
which are similar to our result for triangular proton conhgurations. The distributions of 
eccentricities within the ’hot spot’ model [ 6 ] seem to be centered more towards higher 
values of 62 than in our result. The probable reason is that the ’hot spot’ model lacks any 
constraint on the relative positions of quarks. It is hard to make any direct comparison 
with the other mentioned calculations, however all of them [5-13] predict measurable 
anisotropies in pp collisions. 

6 Summary 

In view of the great interest in understanding proton structure, it is pointed out that 
even simple model ideas concerning distribution of quarks and gluons in protons lead to 
non-trivial multiplicity dependence of initial anisotropies in pp collisions. We illustrate 
this statement with results obtained for phenomenological models intuitively inspired by 
RGPEP in QGD and similar ideas. Expecting new data on high-energy pp collisions from 
the y/s = 14 TeV runs of the Large Hadron Gollider, we believe that detailed studies of 
the role of proton structure models in description of such collisions may shed new light 
on the internal structure of protons. 
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Figure 1: Constant proton density surface for I (quark-diquark) model 
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Figure 2: Constant proton density surface for Y model 
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Figure 3: Distribution of event eccentricities for II, YY and Gaussian-fluctuating type 
(G-fG-f) collisions. The continuous line is reconstructed from the values in particular 
eccentricity bins and represents the normalized probability density of the occurence of 
an event with certain eccentricity. Errors (square roots of variances) follow from basic 
statistics assuming independent events (applies to all other figures). 



Figure 4: Distribution of event triangularities for II, YY and Gaussian-fluctuating type 
(G-fG-f) collisions 
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Figure 5: RMS eccentricities in collision multiplicity bins for II, YY and Gaussian 
fluctuating type (G-f G-f) collisions 



Figure 6: RMS triangularities in collision multiplicity bins for II, YY and Gaussian 
fluctuating type (G-f G-f) collisions 
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Figure 7: Event multiplicity distributions for II, YY and Gaussian-fluctuating type col¬ 
lisions compared to the experimental charged hadron distribution in ^/s = 7TeV pp 
collisions [21] for the pseudorapidity range \ri\ < 2.4 
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